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ABSTRACT 

We describe state-of-the-art MARCS-code model atmospheres generated for a group of A dwarf, G 
dwarf, and late-G to mid-K giant standard stars, selected to photometrically calibrate the Spitzer- 
IRS 4 , and compare the synthetic spectra to observations of HR 6688, HR 6705, and HR 7891. The 
general calibration processes and uncertainties are briefly described, and the differences between var- 
ious templated composite spectra of the standards are addressed. In particular, a contrast between 
up-to-date model atmospheres and previously published composite and synthetic spectra is illustrated 
for wavelength ranges around 8 /im (where the SiO Av = 1 band occurs for the cooler standards) and 
A > 20/im, where the use of the Engelke function will lead to increasingly large discrepancies due to 
the neglect of gravity in cool stars. At this point, radiometric requirements are being met, absolute 
flux calibration uncertainties (1-er) are ~ 20% in the SH and LH, and ~ 15% in the SL and LL 
data, and order-to-order flux uncertainties are ~ 10 % or less. Iteration between the MARCS model 
atmosphere inputs and the data processing will improve the S/N ratios and calibration accuracies. 
Subject headings: instrumentation: spectrographs — methods: data analysis — stars: atmospheres 
— infrared: stars 



1. INTRODUCTION 

The scientific interpretation and modeling of the 
spectra produced by the Infrared Spectrograph (IRS; 
iHouck et alJl2004fl onboard the Spitzer Space Telescope 
IjWerner et al.l2004[) require accurate spectrophotometric 
calibrations, which depend stars with well known envi- 
ronmental and atmospheric properties. We have a unique 
opportunity to use the state-of-t he-art MARCS-code, de- 
veloped by the Uppsala group (Gustafss on et al.1 11975ft 
and modified substantially in the numerical methods and 
input line and continuous opacities, to compute syn- 
thetic spectra of the IRS standard stars. The iteration 
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between stellar models and instrum ent calib rations has 
been described bv lDecin et"aTl ((2000, 2003330), with ap- 
plications to calibrations of the Short Wavelength Spec- 
trometer (SWS) onboard the Infrared Space Observatory 
(ISO). 

In this paper we describe the synthetic spectra based 
on the MARCS-code model atmospheres tailored to a set 
of standard stars that are relied upon for the photometric 
calibration of the IRS. We overview the calibration strat- 
egy, and summarize the primary IRS stellar standards. 
The latest generation model spectra of these stars are dis- 
cussed, presented in comparison to IRS spectra of three 
standards. We finally contrast our synthetic spectrum 
for 7 Dra (K5 III) with widely available Kurucz model 
and composite spectra, and summarize the current pho- 
tometric uncertainties. 

2. THE IRS SPECTROPHOTOMETRIC CALIBRATION 
SCHEME 

The basic strategy for spectrally and photometri- 
cally calibrat i ng th e IRS modules has been outlined by 
iMorris et al.l l|2003h . Each spectral order in the Short 
Low and Long Low (SL and LL, respectively; A/AA ~ 
70 — 140) long slit spectrographs, and the Short High 
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and Long High (SH and LH, respectively; A/AA ~ 700) 
echelle spectrographs is individually calibrated 

2.1. Spectral Calibrations 

Spectral calibrations are determined by observations 
of a combination of smooth, extended celestial sources, 
and emission line sources. The zodiacal light at its maxi- 
mum intensity and extended nebulae pointed off the cen- 
tral source are useful for determining the order widths, 
which reflect the widths of the slits, spectrally imaged 
into the arrays. The order widths are measured at the 
45% crossing point of the illumination profile for each 
order, and are very nearly 5 pixels for each of the SH 
and LH echelle orders, and 36 pixels for the low resolu- 
tion orders. The output signal from the zodiacal light 
measurements can also be used to characterize the pixel- 
to-pixel response variations along the cross-dispersed di- 
rection of the SL and LL orders, but is not sufficiently 
strong for accurately determining these response varia- 
tions in the high resolution orders. Emission line sources 
such as Be stars P Cygni (B2pe I) and 7 Cas (B0.5e IV), 
planetary nebulae NGC7027, NGC6543, SMP083, and 
SMP031, and Saturnian moon Titan all provide strong, 
unresolved lines to determine the wavelength solutions, 
spectral resolutions, and instrument profiles. Moving the 
sources to different locations in the slits allows us to fur- 
ther characterize impacts on the wavelength calibration 
by pointing offsets, and by pixel undersampling of the 
point spread function (PSF). All four arrays undersam- 
ple the PSF, due to the excellent focus of the telescope 
and of the f/12 light beamed onto the IRS arrays, which 
maximizes the sensitivities. However, the undersampling 
leads to various spectral and photometric complications, 
such that (for example) line centroiding accuracies are 
reliable to ~l/5 a resolution element, about half of the 
accuracy achievable for a critically sampled PSF. Mea- 
surements from flight observations of celestial sources su- 
persede the laboratory measurements, carried out in pre- 
launch tests to verify design requirements. 

2.2. Photometric Calibratoins 

Relative spectral response calibrations are determined 
by placing photometric standards at as many locations 
in the slits as is feasible (executed as fine spectral maps), 
processing and combining the data collection events into 
a single 2-D image plane, and removing signature of the 
celestial source by means of the synthetic spectra. Ide- 
ally the spectral flatfield for each array is determined by a 
weighted mean of flats individually derived from different 
stars; initially they are determined for each array from 
a single star. The resultant flatfield is then applied to 
observations of other standards, and the spectral extrac- 
tions are analyzed to determine the flux conversion co- 
efficients (electrons/second to Janskys) and polynomial 
"tuning" coefficients to correct for the effects of diffrac- 
tion losses on point source observations, and systematic 
low frequency residual errors from the flatfielding. The 
flatfields are stored as FITS image files to be used in the 
IRS Basic Calibrated Data (BCD) Science pipeline, and 
the flux conversion and tuning coefficients are stored in 
ASCII-format tables for application to 1-D spectral ex- 
tractions in the post-BCD Science pipeline. 

The implementation of the above scheme requires 
staged efforts. In the Science Verification (SV) phase 



TABLE 1 

IRS PHOTOMETRIC CALIBRATION STARS, AND DERIVED T e g AND 

log g. Uncertainties are in parentheses (cf. Eq. (18) in 
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Source Sp. 


Cal. a 


Range b 


T c ff 


log g 


Type 






[K] 


[cm/s 2 ] 


HR 6688 K2 III 


P 


SH, LH 


4465 (50) c 


2.17 (0.19)* 


HR 7310 G9 III 


P 


[SH],[LH] 


4830 (50) c 


2.58 (0.06) f 


HR 2194 A0 V 


P 


SL1.SL2 


10325 (240) d 


4.09 (0.08) 


HR 7341 Kl III 


P 


LL2 


4570 (50) c 


2.46 (0.08) 


HR 7950 A1.5 V 


S 


SL1,SL2 


9060 (120) d 


3.51 (0.08) f 


HR 7891 A0 V 


s 


SL1,SL2 


10170 (120) d 


4.07 (0.07) f 


HR 6705 K5 III 


s 


SH,LH 


3980 (50) c 


1.06 (0.09) f 


HR 6606 G9 III 


s 


all 


4975 (50) c 


2.90 (0.07) 


HR 2491 Al V 


T 


SH,LH 


10240 (120) d 


4.39 (0.07) f 


HD 105 GO V 


T 


all 


5930 (70) c 


4.31 (0.11) g 



a Calibration types: P = Principal, S = Secondary, T = Testcase. 
Principal standards are observed to derive flatfields and absolute 
flux calibration; secondary standards, which may have limited vis- 
ibilities or less certain stellar parameters, are observed to verify or 
further refine the calibrations. 

b Entries refer to the IRS modules, as follows: SL1 7.5-14.5 (im; 
SL2 5.3-7.5 ^m; LL1 19.5-38.0 /im; LL2 14.0-21.3 fim; SH 9.9-19.6 
/im; LH 18.7-37.2 fim. Brackets denote where calibration stars may 
be used only over certain ranges of the module. 

c Derived from (V — K). 

d Derived from (V — I). 

c Derived from (B — V)-tcmpcraturc relation of Flower (1996). 

f Using [Fe/H] from lCavrel de Strobel et all 119971) . 

s Using [Fe/H] f rom the uvby-[3 — metallicity relation of 
ISchuster fc Nissenl tl989l) . 



only sparse spectral maps of the calibration stars could 
be acquired, and depended on zodiacal light observations 
to spatially fill in the orders. 3 Due to the relatively low 
zodiacal fluxes (in SH, LH, and SL-2nd order), the effects 
of PSF undersampling, and space weather on the detec- 
tors, exhaustive observations of the standard stars con- 
tinue to be carried out in order to meet the radiometric 
requirement of 5 %. The details of the stellar atmosphere 
models are especially important at this stage. 

3. STANDARD STAR SELECTION 

Table ^ summarizes the stars which have been se- 
lected and observed during SV and in IRS Science cam- 
paigns for the purposes of photometrically calibrating 
IRS spectroscopy. The stars were chosen in the pre- 
launch preparatory phase, to meet (as closely as pos- 
sible) specific criteria on the availability of observa- 
tional data to make reasonable estimates of the stel- 
lar parameters (described in the next section), their 
environments, and absence of chromosphcric activity, 
circumstellar dust shell or disk, multiplicity, or pecu- 
liar spectral activity. The IRS operating at low res- 
olution is ~300 times more sensitive than SWS was 
at 10 /im, and consequently the fainter standard stars 
may meet the aforementioned requirements to a lesser 
extent than stars in the ISO calibration programs. 
We make use of observation s from the ISO ground- 
based preparatory program IjHammerslev et aLl l!99St 

3 The zodiacal light itself is unreliable in the spectral dimension 
due to the likely presence of solid state features and latitudinal 



varia tions in dust temperature and grain properties IReach et alJ 
2003). 
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iHammerslev fc Jourdain de Muizoia |2003|) and detailed 
theoretical work where possible. In particular, HR 6688, 
HR 6705, HR 7310, and HR 2491 in Table Q] have 
been previously modeled bv lDeciiiet al. (2003a|g) using 
ground-based and 2.4 - 12fim ISO-SWS spectroscopy. 
Generally, G-K giants are used to calibrate the LH, LL, 
and SH modules, and the A dwarfs to calibrate SL. 
This balance is struck by the higher potential for debris 
disks around the A stars, and non-photospheric molec- 
ular and dust layers around late-type M giants. The 
selection of the G and K standard stars was nevertheless 
made on a very restricted basis: several infrared stud- 
ies based on the CO Av = 2 lines (at 2.3 /im) and the 
CO Ad = 1 lines (at 4.6 /im) have revealed a thermal di- 
choto mic structure in the out er layers of later-type stars 
(e.g. Wiedemann ct al. 1994). A two-component struc- 
ture — consisting of the traditional chromosphere and a 
radiative equilibrium region mediated by molecules with 
CO cooling dropping the temperature — with physically 
distinct areas of hot and cool material at the same al- 
titude is thought to better represent the cool stellar at- 
mosphere than do existing homogeneous one-component 
models. For our purposes, we have selected stars which 
are thought to belong to the group of "quiet" stars for 
which the radiatively-cooled regions very largely dom- 
inate the stellar surface, so that the single component 
RE atmospheric models are most likely to be valid. 

4. REFERENCE SEDS — MARCS SYNTHETIC SPECTRA 

A new set of theoretical reference SEDs has been cal- 
culated using the MARCS and turbo-s pectrum code 
IjGustafsson et alJ ll975: Ple z et aLlll992l and further up- 
dates), using the same physical input parameters of 
atomic and molecular equilibrium constants, solar abun- 
dances, continuous opacity sources, etc., as the ones de- 
scribed by Decin (2000). For line opacities in the IRS 
spectral range, a database of infrared lines with atomic 
and molecular transitions (CO, SiO, H 2 0, OH, NH, HF, 
HC1, CH, and NO) has been prepared. Referen c es an d 
discussion of the inputs can be found in iDecinl (2000) . 
The standard assumptions of homogeneous stationary 
layers, hydrostatic equilibrium and LTE were made. En- 
ergy conservation was required for radiative and convec- 
tive flux, where the energy transport due to convection 
was treated through a local mixing-length theory. The 
turbulent pressure was neglected. The Rosseland optical 
depth scale has been chosen to span values from 6xl0 -8 
to 300 in order to ensure the fulfillment of the diffusion 
approximation adopted as lower boundary condition for 
all frequencies, and to minimize the the number of fre- 
quencies for which the surface layers are still optically 
thick. 

The computed theoretical atmosphere model and syn- 
thetic spectrum depend on a large number of input pa- 
rameters, the main ones being the effective temperature 
T c ff , the gravity g, the microturbulence and the chemical 
composition. In case of a spherical symmetric geometry, 
one also has to provide either the mass M or the radius 
R. the s tellar parameters as determined bv lDecin et alJ 
(2003alcl) for stars in common with the ISO-SWS calibra- 
tors may be adopted. However, since we aim to set up a 
spectral response calibration independently, we prefer to 
make initial estimates of T c s and log g [cm s~ 2 ] from pho- 
tometric colors (see below). Since almost no information 
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FlG. 1.— High resolution spectra of HR 6688 (top) and HR 6705 
(bottom), with the MARCS model SEDs overplotted in red. S/N 
ratio estimates are labelled. The insets show regions of the IRS 
spectrum around OH lines, compared to the model (in red, offset 
by -0.1 Jy). 



is available for the microturbulence £t, a value of 2 km 
s _1 was assumed. Nevertheless, it should be noted that 
since the final calibration of the ISO-SWS relative spec- 
tral response functions (RSRF) was also based on MARCS 
model atmospheres, the IRS spectral response calibration 
and methods to ite rate between the input parameters and 
spectroscopy (e.e. ADecm et alJ l)2003alb| |q)'l are linked to 
the SWS-RSRFs. Moreover, first estimates of C, N, and 
O abundances and the 12 C/ 13 C-ratio — being crucial 
for a correct computation of the molecular opacities in 
K-giants — are ado pted from recent an alyses of the ISO- 
SWS calibrators bv lDecin et alJ l)2003a|) . T eff can be cal- 
culated for K and M giants directly from (V—K), and for 
A-type dwarfs from (V — I), usin g (semi)-e mpirical color- 
temperature relations (e.g. .lBes sell et aLl |1998|) . Such a 
calibration is e.g. given by[Bessell ct al. (199§) who de- 
rived a polynomial fit between T e g and (V—K) or (V—I) 
from (a) the infrared-flux method (IRFM) for A-K stars, 
and (b) Michelson interferometry for K-M giants. For 
that purpose, the (V-K) and (V-I) Johnson or 2MASS 
colors were converted to th e Johnson-Cousins-Glass sys- 
tem IjBessell fc BrettJll988j) and correcte d for interstellar 
extinction using Ay = 0.8 mag kpc -1 iJBlackwcl l et alJ 
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ITOflfl). and E (V - K ) = Ay/1.1 and Aj = 0A8A V 
ljMathisll990j) , with the distance calculated from the Hip- 
parcos parallax n. Whenever the de rived color est imates 
were outside the ranges specified by Bcssell et al. (1998) 
in their Table s 7-8, the co lor-temperature relation as 
determined bv iFlowerl l)1996ft was used. 

In order to estimate the gravity, one needs the radius R 
and the stellar mass M. The first parameter is assessed 
from (K, BC^, tt, T off ((V - K) )) for K giants or (V, 
BCy, n, T cff ((V - J) )) for A-G dw arfs. The BCV bolo - 
metric corrections are derived from lBessell et aLl([ l998') 
whene ver appropriate, otherwise the BCV data of lFlowerl 
( 1996) were used. The uncertainty on the bolometric cor- 
rectio n is assumed to be 0.05. Mbol,© is assumed to be 
4.74 IjBessell et al.lll998j) . Mass values for the stars in 
our sample are estimated from evolutiona ry tracks with 
appro priate metallicity as calculated by IGirardi et al.l 
(2000). The estimated mass depends critically on the 
assumed metallicity, which has been adopted from liter- 
ature references (see Tabled . Where no information was 
available, a solar metallicity was assumed. The resultant 
gravity for each star is listed in Table ^ 

5. IRS SPECTRA OF REPRESENTATIVE STANDARDS 

In this section we compare reduced spectra of repre- 
sentative standards to their synthetic SEDs, to assess 
the models and calibrations. The observations have been 
processed to BCD products using the SSC pipeline ver- 
sion S9.5.0, and were response calibrated with the lat- 
est flatfields, except for HR 6688, which was calibrated 
from observations of Sirius and Vega. Prior to extrac- 
tion of low resolution spectra of HR 7891, the flatfielded 
BCDs were treated for removal of background sky emis- 
sion. Spectra were extracted and flux calibrated with the 
SSC pipeline, from BCDs with the star near the nod and 
slit center positions, and then spectra were sigma-clipped 
and combined. Photometric uncertainties are discussed 
in Sect. 
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Fig. 2. — Background-corrected low resolution spectrum of 
HR 7891, with the MARCS model spectrum overplotted in red. The 
means and standard deviations of the ratio between observations 
and the model SED in each spectral order are indicated. The 
higher spread of observed LL-1 fluxes about the model is due to 
the 37.5— 40 /^m range where throughput response is very low, and 
where a soft filter cutoff allows some contamination with second 
order light from ~20^tm. 



5.1. Bright standards HR 6688 and HR 6705 

The comparison between the IRS spectra and the 
model of HR 6688 is shown in the upper panel of Fig. Q] 



The data quality is sufficiently high to make plausi- 
ble detections of spectral lines with predicted line-to- 
continuum ratios as low as 1% (see Fig^inset displaying 
OH Av — spectral features) . The shape of the observed 
continuum is in excellent agreement with the model (re- 
lying intrinsically on the accuracy of the Sirius model). 
The high resolution spectrum of HR 6705, the primary 
flux calibrator for the ISO-SWS, is shown in the lower 
panel of Fig. ^ The region between 10— 15/xm is not 
plotted due to saturation over portions of echelle orders 
15—20. The remaining spectrum is calibrated with the 
default calibrations, that is, HR 6688-based flatfields and 
tuning factors. The shape of the continuum is again in 
excellent agreement with the model, and the inset plot 
again demonstrates the potential for detection of weak 
(molecular) spectral features, and improvements to pro- 
cessing methods. 

5.2. HR 7891 

The spectrum of HR 7891 plotted in Fig. |3 spanning 
two orders of magnitude in flux densities, is representa- 
tive of the agreement between independently calibrated 
low resolution spectroscopy of this A dwarf secondary 
stardard star, not previously observed at these wave- 
lengths, and its MARCS model atmosphere. The back- 
ground sky was subtracted with off-pointed subslits prior 
to spectral extraction. The ratio the model and observa- 
tions show a ~10 % error in the absolute flux calibration 
of the SL-1 portion, but otherwise excellent agreement, 
and that the star has no detectable thermal excess from 
a potential debris disk. 

It is helpful to point out the difficulties imposed by 
background emission, increasing towards the ecliptic 
plane to typical levels of several tens of MJy sr _1 . Nor- 
mally this can be corrected with an off-pointing, as done 
for the observations plotted in Fig. [21 If this is not avail- 
able, then a coarse correction can be applied using a 
Spitzer-centric COBE/DIRBE thermal model of the zo- 
diacal dust on the date of the observation, using appro- 
priate scaling factors for the solid angles of the individ- 
ual slits. The uncertainties associated with the zodiacal 
light model (including the likely presence of silicate dust 
bands), and with the aperture solid angles, limit the ac- 
curacy of the corrections to ~20% at any wavelength. A 
residual tilt may also remain from the fact that both ex- 
tended and point source signal may fall through the slits, 
while calibrations are based solely on point sources (such 
that corrections for diffraction losses would not be ade- 
quately compensated by simple scalings using a zodiacal 
light model). 

6. COMPARISON WITH TEMPLATED SPECTRA 

The photometric calibration schemes of the ISO and 
Spitzer instruments have been partly based on templated 
datasets, such as those constructed by M. Cohen and col- 
laborators (e.g.. lCohen et afll2003(K with both observa- 
tional and theoretical components to them. The Cohen 
datasets can be divided into three categories: (a) a col- 
lection of absolutely calibrated photometric and spectro- 
scopic data, merged to create a continuous spectrum with 
uncertainties traceable to a specific group of stars such 
as a Tau and Vega; (b) a Kurucz model adopted at a lit- 
erature concensus of T c ff , log(g), and metallicity, then fit 
to available photometry; and (c) a template constructed 
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Fig. 3. — Comparison between reference SEDs of HR 6705. The composite from Coh en et al.l 11993) is based on different photometric 
data (orange dots) and spectral fragments spliced together. The used data sets are indicated by the solid black line, and quoted uncertainty 
limits by the dashed lines. The input parameters for the MARCS (red) and Kurucz (blue) synthetic spectrum are indicated. In the bottom 
panel, the wavelength range observed by the IRS is enlarged, and the spectral coverage of the SL, LL, SH, and LH are indicated. 



from a star from (a) and/or (b) with the same spectral 
type, and scaled to the observed photometric data. Stars 
with the same spectral type may exhibit a different abun- 
dance pattern, affecting e.g. the SiO band strengths and 
continuum levels. Products in the (a) group when avail- 
able are helpful to the IRS since they rely on no theoret- 
ical assumptions, except for the use of an Engelke func- 
tion at A > 20 /im, and some products in the (c) group 
are also useful for "testcase" standards which can be ob- 
served to monitor the external photometric stability of 
the instrument in conjunction with pointing performance 
of the spacecraft, AOT repeatibility, and statistical mea- 
sures. A number of stars not listed in Table n are being 
observed in normal operations for these activities, and 
we rely on products in groups (a) and (c) for stars whose 
individual value is their membership in the ensemble. 

For the stars listed in Table ^ we re ly explicitly on 
the MARCS synthetic spectra for the detailed calibration 
analyses, noti ng that while a composite is available for 
HR 6705 (Cohen et al. 1996), it makes use of an En- 
gelke function extrapolation beyond 18.5 /zm, which we 
prefer to avoid (see next paragraph). For a wider (sec- 
ondary) set of stars, which are generally too uncertain 



in their stellar or environmental parameters to initially 
justify the computational resources of MARCS modeling, 
we can use composite observations templated by spectral 
type and luminosity class, or Kurucz synthetic spectra 
with the understanding that they are calculated assum- 
ing a plane-parallel geometry, solar abundances (scaled 
to the metallicity), and a microturbulent velocity of 2 km 
s _1 , while these are all free input parameters for the 
MARCS atmosphere models. Note also that some of the 
line lists for diatomic molecules used in Kurucz models 
are fifteen years old, incomplete, and do not reproduce 
line positions very well for high v and high J. The stars of 
this group are typically used for trending and checks of 
order-to-order and module-to-module calibrations, and 
response to high and low flux point sources. 

To illustrate these differences over the mid-IR spec- 
tral range, we show three datasets created for HR 6705 
in Fig . |21 The composite constructed by iCohen et alJ 
(1996) is plotted in black, with the photometric points 
and the various spectral fragments indicated. Of great 
importance to the IRS occurs at A > 18.5 /xm, where 
fluxes are ap proximated with the Engelke function 
l)Engelkei ll992). which offers an analytical implicit seal- 
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ing of a semi-empirical solar atmospheric temperature 
profile to differing effective temperatures. This analyti- 
cal (plane-parallel) approximation should be valid for the 
2-60 /im continuum for giants and dwarfs with effective 
temperature 3500 K < T c ff < 6000 K, where scattering of 
electrons from H - dominates the continuum op acity. By 
fitting this function to a set of standard stars, lEngelke 
(1992) concluded that the estimated probable error in 
estimated flux is ±3 % below 10 /im, up to ±5 % in the 
vicinity of 25 iim, and 6 % at 60 (im. The main spec- 
tral difference between the MARCS and Kurucz theoreti- 
cal spectra is the absence of the SiO fundamental (around 
8/im) and first-overtone (around 4/im) lines in the Ku- 
rucz spectrum. Smaller differences do occur due to the 
use of a different abundance pattern. With an exten- 
sion d = (R(t ioss = 10~ 5 )/i?(T ross = 1)) - 1 being only 
3.7%, the spectral differences between a plane-parallel 
and spherical geometry are only marginal. 

The lowermost plot in Fig.|3focuses on the wavelength 
ranges of the SL, LL, SH, and LH modules. Two main 
differences between the 3 data-sets are of major concern 
for IRS calibrations: (1) around 8/im where the SiO 
Av = 1 lines occur, and (2) for A > 18.5 /jm. 

(1) Comparison to the LRS data indicates that the 
abundance pattern used as input for the MARCS model 
and synthetic spectrum calculation should be altered, in 
particular the oxygen (and carbon) abundance. With 
the moderate resolution of IRS (~ 600) and the SL mode 
only starting around 5 /xm (being in the middle of the CO 
Av = 1 band for cool giants), it will be very difficult to 
constrain the stellar parameters (and in this case more 
specifically the oxygen abundance) from the IRS data. 
This situation can be avoided by relying on A dwarfs 
and G giants in this wavelength range. 

(2) At LL and LH wavelengths, we see a clear differ- 
ence between the slope of the MARCS, Kurucz, and Cohen 
(Engelke) spectrum. With the Engelke function being 
a plane-parallel approximation, we have compared this 
function with MARCS plane-parallel spectra for T e ff be- 
tween 3500 and 6000 K for different values of the gravity. 
The best agreement occurs for T c g- — 6000 K, which is not 
surprising since Engelke's function is based on a scaled 
solar model. For lower effective temperatures, the role of 
the gravity becomes more important: a higher gravity re- 

4 The low background limit is important even for the high resolu- 
tion modules. Automatic background correction is not performed 



suits in a higher opacity (from the free-free transitions of 
H _ ), a more compact object and a smaller temperature 
gradient due to the requirement of flux conservation. As 
a consequence, the infrared continuum slope (A > 2 /im) 
is less steep for a higher gravity model. From a com- 
parison with the MARCS cont inua, we may co nclude that 
the uncertainties quoted by lEngelkd l|1992j) do not ac- 
count well for the influence of gravity. For HR 6705, the 
difference in slopes between the MARCS model and the 
composite spectrum (using an Engelke function at A > 
18.75 /an) produces systematic underpredictions of fluxes 
in the composite by 4.2% at 30 /zm, and 6.3% at 70 pm. 
These differences are easily within grasp of the sensitive 
IRS and the MIPS detectors. 



7. SUMMARY OF THE SPECTROPHOTOMETRIC 
UNCERTAINTIES 

First, we estimate that for well-pointed observations 
of point sources, the relative photometric uncertainties 
within any spectral order generally meet the 5% radio- 
metric requirement, over the spectral ranges committed 
to observers. The exception is the two reddest orders 
of LH (34.2 - 37.2 /im), where throughput response to 
HR 6688 is low, and will be improved with HR 6705 
calibration observations by the time that this paper is 
published. It must also be noted that the ranges of 14.1- 
15.4 lim in SL 1st order and 38-42.4 /zm in LL 1st order 
may be contaminated with light from around 7-8 and 
19-21 pm, respectively, due to a weakness in the trans- 
mission cutoff of the filter in the 2nd spectral orders. The 
level of contamination depends on the color of the source 
in the slit. 

The photometric performance of the IRS is very sensi- 
tive to pointing effects, due to the sizes of the slits and 
the PSF widths. Combined with the effects of PSF un- 
dersampling, the absolute flux calibration is estimated at 
this point to be uncertain by 20% in the SH and LH data 
overall, and 15% in SL and LL data. Generally, order- 
to-order calibrations are 10% or less, for point sources 
in the low background (or background-corrected) limit 4 
and well-placed in the center or nod positions of the slits. 
The errors are easily recognized, and can be mitigated by 
scaling to photometry, where available. 

in the SSC pipeline. 
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